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The geometry of impulsive pp-waves is explored via the analysis of the 
geodesic and geodesic deviation equation using the distributional form 
of the metric. The geodesic equation involves formally ill-defined prod- 
ucts of distributions due to the nonlinearity of the equations and the 
presence of the Dirac (^-distribution in the space time metric. Thus, 
strictly speaking, it cannot be treated within Schwartz's linear theory 
of distributions. To cope with this difficulty we proceed by first regu- 
larizing the 5-singularity, then solving the regularized equation within 
classical smooth functions and, finally, obtaining a distributional limit 
as solution to the original problem. Furthermore it is shown that this 
limit is independent of the regularization without requiring any addi- 
tional condition, thereby confirming earlier results in a mathematical 
rigorous fashion. We also treat the Jacobi equation which, despite being 
linear in the deviation vector field, involves even more delicate singu- 
lar expressions, like the "square" of the Dirac ^-distribution. Again 
the same regularization procedure provides us with a perfectly well 
behaved smooth regularization and a regularization-independent dis- 
tributional limit. Hence it is concluded that the geometry of impul- 
sive pp-waves can be described consistently using distributions as long 
as careful regularization procedures are used to handle the ill-defined 
products. PACS-numbers: 04.20. Cv, 04.20.-q, 02.20.Hq, 04.30.-w 



I. INTRODUCTION 



Plane fronted gravitational waves with parallel rays (pp-waves) are spacetimes 
characterized by the existence of a covariantly constant null vector field, which 
can be used to write the metric tensor in the formEl 

ds 2 = H(u,x,y)du 2 — dudv + dx 2 + dy 2 , (1) 

where u, v is a pair of null coordinates (u = t — z, v = t + z) and x, y are transverse 
(Cartesian) coordinates. In this paper we shall deal especially with impulsive (in 
the diction of Penroseo) pp-waves where the profile function H is proportional to 
a (5-distribution, i.e., takes the form H(u,x,y) = f(x.y)S(u) 1 where we leave the 
(smooth) function / of the transverse coordinates arbitrary for the moment. This 
metric is flat everywhere exept on the null hypersurface u — 0, where it has a 
(5-shaped "shock". 

Such spacetimes arise most prominently as ultrarelativistic limits of black holp ge- 
ometries as first derived by Aichelburg and Sexl for the Schwarzschild casetJTJ. On 
the other hand Penrose has given a more intrinsic description of such spacetimes by 
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his "scissor and paste" approach]], which essentially consists of glueing together two 
pieces of Minkowski spacetime along the null hypersurface u = with a shift in the 
v-direction. A similar idea was used by Dray and t'HooftEl, who introduced a coordi- 
nate shift along geodesies in Minkowski space time to rederive the Aichelburg-Sexl 
geometry as well as Penrose's junction conditions from the field equations. 
The philosophy of the present work is somehow complementary. We take the 5- 
shaped metric literally and try to explore the properties of this geometry via inves- 
tigation of geodesies and the geodesic deviation, thereby following the approaches 
of Ferrari, Pendenza and Veneziandj and BalasinEI. 

The main purpose of this work is to deal with the singular, i.e., distributional quan- 
tities in a mathematically rigorous fashion. A thorough analysis of the geodesic 
equation shows that it involves ill-defined products of the (5-distribution with the 
step function. This difficulty can be circumvented by a proper regularization pro- 
cedure providing us with a perfectly well behaved smooths-approximation, which 
has a distributional limit coinciding with the earlier resultsa-EI. However, we neither 
have to impose "multiplication rules" like 59 = (1/2)5, nor to use any additional 
requirements as the constancy of the norm of the geodesic 's tangent vector across 
the shock. In fact, this property comes out as a result in our approach. The details 
of this calculation are given in Sec. |HJ 

Further investigations show that the Jacobi equation involves ill-defined terms of 
an even worse type, such as the "square" of the Dirac (^-distribution and the square 
of the step function times the ^-distribution. However, our regularization strategy 
again provides us with a smooth approximation and a well behaved distributional 
limit, even in this case, where it seems to be hopeless to use "ad-hoc extensions" 
of Schwartz's linear distribution theory such as special "multiplication rules". We 



give the details of the rather lengthy calculation in Sec. HI. Our distributional 
"solution" of the Jacobi equation fits perfectly well in the (hcuristically) expected 
picture, showing the consistency of our approach. Finally, we make some comments 



and give an outlook to future work in Sec. [V 



II. GEODESIC EQUATION 



We start with an impulsive pp-wave metric of the form 

ds 2 = f(x l ) 5(u) du 2 -dudv + (dx 1 ) 2 , (2) 

where x % (i ~ 1, 2) denote the transverse coordinates. It is straightforward to derive 
the geodesic equation. The nonvanishing Christoffel symbols are 

Ku = -fs, ri u = ~ dj 5, I™ = = -dj 5 , (3) 

where we have denoted the partial derivatives of / by dif and the derivative of the 
(^-distribution by 5. Hence we get the equations 

u" = , 

v" = f6 + 2(d i fx i ') u'5, 
x^' ^^fu' 2 5, (4) 

where ' denotes the derivative with respect to an affine parameter and summation 
over i is understood. We use the first equation to introduce u as a new affine 
parameter (there excluding trivial geodesies parallel to the shock hypersurface) to 
get 
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v(u) = f(x*(u))5(u) + 2 (difix^u)) x l (u)) 6(u) , 
x i (u) = ±d i f(x i (u))8(u) , (5) 

where ' again denotes the derivative with respect to u and we have inserted all the 
dependences explicitly. Equations (||) form a system of three coupled, nonlinear 
ODEs of second order in the vectorspace V of distributions. For u ^ all the right 
hand sides vanish, which is clear from the form of the metric tensor, and we expect 
the geodesies to be broken, possibly refracted straight lines. However, if we take a 
closer look at system (|J) we see immediately that the first equation cannot be taken 
literally in the sense of distributions as the terms x % 5 involve ill-defined products 
of the ^-distribution with the step function. 

To analyze the situation in some more detail we integrate the last two equations 
using the (distributional) identity f(x l (u)) 6(u) = f(x l (0)) S(u) to get 

x l (u) = initial values + - d% f(x l (0)) u+ , (6) 

where we have denoted the "kink" -function uO(u) by u+. Note, however, that dis- 
tributions can only be multiplied by C°°-functions, whereas it's not clear a priori 
that the solution x l (u) will be smooth; in fact as suggested by equation @ and as 
shown later in our calculations the solution will not even be differentiable at u = 0. 
If we still try to go on by brute force it comes even worse: inserting (Q) into the first 
equation (||) we see that the term x % {u) 5(u) gives-rise to the ill-defined product 95. 
(For some further comments on this product see£3, p 21.) 

To overcome the undefincdness described in detail above, we apply a careful reg- 
ularization procedure. More precisely we regularize the (5-distribution by a (stan- 
dard) mollifier or model <5-net, i.e., a net p ei defined as follows. Let p be a smooth 
function with support contained in the interval [—1,1] and J p = 1; now put 
p e (x) — (1/e) p(x/e), (e > 0). Hence system (^) takes the regularized form 

v e (u) = f(xi(u))p e (u) + 2dif(xi(u)) ±\{u) Pe {u) , 

K{u) = \d l f{<{u)) Pe {u) ■ (7) 

Next we choose initial conditions in u = —1 (i.e. "long before the shock"), more 
precisely 

v e (-l) = v 0l xl(-l) = 4 , 

«e(-l) = vo, = 4 . (8) 

for all e. Note that choosing initial conditions in u = would mean to "start at the 
shock" and one cannot expect to (and indeed does) end up with a regularization 
independent result in this case. 

Due to the regularization procedure we now have to deal with the fully nonlinear 



character of the last equation of (|7]). However, as shown in appendix IV the spe- 
cial form of the right hand side of the equation guarantees global existence of the 
solutions for small e. These are given implicitly by 

v e (u) = v + v (1 + u) + {0 * 9 * v e ) (u) , 

x\(u) = 4 +4 (1 + u) + (d * 6 * 4) (u) , (9) 
where "*" denotes convolution. 
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We are going to calculate the (distributional) limits of the solutions (^) as e tends 
to zero. Since distributions supported in an acute cone form a convolution algebra 
(where, in particular, convolution is a separately countinuous operation) it suffices 
to calculate the limits of the right hand sides of (Q); the distributional limits of the 
solutions (H) are then computed simply by integration. 

We begin with the latter two equations of system (fjj) and choose a test function ip. 
We have to calculate the limit of 

1 

If I f 

dif(xl(u))p e (u)ip(u)du = - dif{xl(eu))p(u)ip(eu)du. (10) 



2 / J y ty " r v ,rK ' 2 



Since (for small e) x e (u) is bounded uniformly on compact sets (see Appendix IV) 
we can use Lebesgue's dominated convergence. Hence the only term we have to 
compute is 

eu v 

^■ V[ ^ X e{ eU ) = um ( X + x ll 

— e — e 

eu v/e 

= 4+4+^lim J J dif{x\(es)) p{s)dsdv . (11) 

-e -1 



l//l< 11^/11 IIpII IJ dsdv 

For e — y the last integral gives zero, since in the limit the range of integration 
only covers a set of zero measure. Hence we get (as expected) 

= \di f(xl(u))p e (u) - \ d t f(xl + ±1) 5(u) , (12) 
within distributions. 

Next we turn to the "critical" first equation in (0) and again calculate the limit of 
the right hand side, i.e., of the expression 

v e (u)=f(xi(u))p £ (u) + 2d i f(xi(u))xl(u) p e (u) 

= (f{x\{u)) Pe {u))- + d % f{x\{u))x\{u)p t {u) . (13) 

The limit of the first term is easily seen to be f(x l + x z ) 6(u), and we are only left 
with the second term which contains the "critical products" . Inserting (0) and (0) 
respectively, we have 

d t f(xi) x\ Pe = f(4) {±1+8* x\) Pt 

= d % f(xl) ±1 Pe + d t f{ X \) (9 * [1 8i /(**) P e] ) Pe ■ (U) 

The limit of the first term again is easily seen to give x l Q di f(x l + x l ) S(u), and we 
are finally left with the task of computing the limit of the convolution term, i.e., 
the expression "#<5". Using 

u u/e 

(B*[8 i f(x i e )p e ])(u) = ( dif(xi(s))p e (s)ds = f dif(xi(es))p(s)ds (15) 
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and again denoting by tp a test function we get 

1 u 



Jim(^ / d t f{xl{eu)) [ j di f(x\(es)) p(s) ds] p(u) <p(eu) du 



l(d t f(xl + il))\(0) I [ I p(s)ds}p(u)du 



2 



1/2 



= \(d i f(xi+x i ))%(p). (16) 
Collecting things together we have 

iim« e = /(4 + 4)j + 49 i /(4 + 4)«5 + i(a i /(4 + 4)) 2 s (n) 

e— >0 4 

within distributions. 

At the end of this section let us give a summary on what we have done so far. We 
have regularized the distributionally ill-defined geodesic equations (|^) by replacing 
the ^-distribution by a (generic) class of mollifiers to obtain ([?]). These equations, 
although nonlinear, provide us with global solutions v e and x\ for small e (as shown 
in appendix [V), implicitly given by (^). Using the latter formula we have shown 



that the smooth solutions have a regularization independent distributional limit 
given by 

v e (u) -^v Q + v (l + u) + /(0) 6{u) + d t /(0) (4 + \di /(0) ) u+ , 
xi(u)^xl+xi(l + u) + ^d i f(0)u + , (18) 

where we have used the abbreviation /(0) = f(x t + Xq). 

Hence viewed distributionally the geodesies are given by refracted, broken, straight 
lines as expected. Of course equations (|l8|) coincide with the earlier resultsErEl How- 
ever from the point of view of our approach the (deeper) reason why here the "rule" 
95 = (1/2) 8 used bycl (which in fact coincides with the "determination of the point 
value" 8(0) — 1/2 used bya) leads to a physically reasonable result is the following: 
The geodesic equations involve only one singular object and hence the <5's as well 
as the 0's appearing above share the same root: namely the i5-shaped wave profile. 
Hence, when regularizing the equations (which in fact corresponds to the physical 
idea of viewing the impulsive wave as an idealized sandwich wave) both factors of 
the ill-defined product naturally involve the same regularization which immediately 
leads to the (regularization-independent) result p t J p t — > (1/2)5. 
We thus conclude that the geodesic equation can be treated consistently by reg- 
ularization, leading to a regularization-independent distributional result. This, of 
course, is only possible due to the relatively mild character of the singular terms 
which allows for a distributional limit of the solutions to the regularized problem 
at all. However, we shall see in a moment that even in the considerably more com- 
plicated case of the Jacobi equation our strategy can be applied successfully. 
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III. JACOBI EQUATION 



In this section we solve the Jacobi equation for an impulsive pp-wave. To keep 
formulas more transparent we restrict ourselves to the axisymmetric case. More 
precisely we restrict the function / of the transverse coordinates x l = (x,y) in the 
metric tensor (pi) to depend on the two-radius r — \/x 2 +y 2 only and work entirely 
within the y = O-hypersurface (initial conditions Xq = ya — 
more we take initial values Vq 
simplify to 



T 2 



yo)- Further- 

= xq. With these assumptions equations ( |l8| ) 



v e (u) 

x\(u) = x € (u) 
x 2 e (u) = y e (u) 



r„(l • +- A f'{x Q fu\ , 



X Q + 2 f'( x o) u + i 



= , 



(19) 



where /' denotes the derivative with respect to the single variable r = x. Hence we 
shall deal with a geodesic tangent vector of the form 



/I 



T a {u) 



vo + f(x )S(u) 

lf(x )e(u) 



\f'{x ) 2 6{u) 



\ 



(20) 



f(x )6(u) + (l/4)f'(x ) 2 9(u) and 



where we are going to use the abbreviations A 
B := (1/2) f'(xo) 8(u) for its components. 

Our next task is to compute the explicit form of the Jacobi equation D d ^ = 
-Rl cd T b T d N c for a vector field N a {u) = (N u (u),N v (u),N x (u),Ny(u)) over the 
geodesic. After some (tedious) calculations we end up with the following form of 
the system 

N u = , 

N v = 2[N x f'S] ' - N x fS + [N u f5]"- N u f"B 2 5 - N u f'BS , 



N x = [tfUfl + 1 N xf,f} g + 1 f, N ug 

N y = , 



(21) 



where we have suppressed the dependence on the parameter u and the variable a;. 
Equations ( pl| ) form a system of four coupled ODE's linear in the components of 
the vector field iV a but nonlinear in the derivatives of the metric. From the fact 
that B involves the step function we see immediately that (in the second equation) 
we again have to deal with distributionally ill-defined expressions, but now of even 
worse type than before. Indeed the term BS is proportional to the "square" of 
the Dirac (5-distribution, and the term B 2 5 involves an expression u 9 2 8". Note, 
however, that the critical terms arise from the second covariant derivative, where 
some of the Christofiel symbols get multiplied, and not from the Riemann tensor 
which components are just proportional to the ^-distribution. To overcome these 
problems we apply the same regularization procedure as in the case of the geodesic 
equation. In particular, we use the regularized geodesic tangent vector 



/I 



V o + * V e 



(22) 



Vo 
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where 

v e (u) = f(x e (u)) p e (u) + 2f'(x e (u)) x e (u) p e (u), x e (u) = (1/2) f'(x e (u)) p e (u) and 
from now on we use the new abbreviations A e := 9 * v e and B e := 9 * i e . Denoting 
the regularized Jacobi field by N?{u) = (N?(u), N"(u), N* (u) , N* (u)) , system (§!]) 
takes the regularized form 

= , 

ftv = 2 [N?f{x e )p e ] ■ - N?f(x e )p\ + [N?f(x e )p £ ]"- 
N?f"(x e )[9 * x e ] 2 Pe - N?f'(x e )x ePt , 

= [N?f'(x e ) + I N*f"(x e )} Pe + \ f'(x e )N?p e , 

= . (23) 

To maintain transparency of formulae we choose appropriate and simple initial 
conditions on the Jacobi field at u — — 1, i.e., 

K(-\) = (o,o,o,o) , 

N?(-l) = ( ai b,0,0) , (24) 

for all e. Note that this corresponds to a focal point at u = —1, and that "nearby" 
geodesies have relative initial velocities only in the u- and w-direction. 
The first and the last equation of system ( p3| ) are easily solved to give 

N^{u) =a{l + u) , 

N*(u) = 0. (25) 



For the remaining, more complicated equations of system ( |23| ) we apply the same 
strategy as in the preceeding section. Since the equations are linear in the com- 
ponents of the deviation vector field we obtain globally defined (smooth) solutions 
which, due to the initial conditions (0), arc implicitly given by 

A^(tt) = b (1 + u) + [6 * 6 * N^] (u) , 

N^(u) = [9*6 * Nf] (u) . (26) 

Again it suffices to compute the distributional limits of the right hand sides of 
equation ( ^3| ) , since by continuity of the convolution we immediately get the limits 
of N? and Nf. 



We start with the third equation of the system (|23|). The main problem we have 
to face here is due to the fact that the unknown function Nf (which, in the limit, 
we cannot even expect to be continuous at u = 0) appears on the right hand side. 
Inserting the initial conditions ( pi} ) and the solutions ( |25| ) we get 

Nf(u) = af'(x e (u))p e (u) + y'(x e (u))a(l + u)p e (u) + l;Nf(u)f''(x e (u))p e {u) . 

" ' ^ . , 2l „ . 

j v v 

lie Hie 

(27) 



The distributional limits of the first two terms in equation fl27j) are easily seen (using 
similar methods as in the previous section) to be 
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lim I e — af (xo) 5 , 

limI7 e = ~ af(x ) (5-5)-]- af'(x )f"(x ) 5 . (28) 

However, the last term in equation ( p7| ) requires a more detailed analysis. Inserting 
again 6*9* Nf for the unknown function Nf and denoting by ip a test function we 
obtain 

III e (u) = 

1 11 S 

~e J p(eu)f"(x f ieu))p(u) J j \af (x 6 (cr)) + \ N?(er)f" \x e (er))] p(r) dr ds du 

-1 -1-1 

1 US 

+ J (p(eu)f"(x e (eu))p(u) J J f'(x e (er)) p(r) dr ds du 
-l -i-i 

1 us 

+ ^ae J ip(eu)f"(x e (eu))p(u) J J f'(x e (er)) rp(r) dr ds du . (29) 

-l -i-i 

The last term in equation ( p9| ) goes to zero, since ip, /', /" are smooth functions and 
x e is bounded uniformly on compact sets (see appendix IV). For the same reasons 
the second summand in equation (^9|) approaches the limit 

ia^(0)/'(x o )/"(x o ), (30) 

whereas the first line also vanishes in the limit, since N'f(er) is bounded on the 
compact range of integration. (Note, however, that lim N'f is not continuous at 0!) 
In some more detail the latter argument can be seen from 

u s 

N*{eu) = J J Lf\x e (er))ep(r) + ^af\x e (er))[p(r)+rep(r)}^ drds 

-l -l 



eu s/e 

1 



N*(er)f"(x e (er))p(r)drds, (31) 

-c -1 

where we have now used equation ( p6[ ) for the third time. Hence we have for u < T, 
where T is some constant, and for fixed e 

eti s/e 

I Nf{eu) | <C + C [ I | Nf(er) \ \ p(r) \ drds 



< C + 2C" j | N?(er) \ \ p{r) \ dr , (32) 

-i 

where C = C(T) and C denote constants. From Gronwall's inequality we gain 

I K(eu) \< Ce 2C ' S -^ p[r) \ dr , (33) 
which is bounded by a constant for u < T. 
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Collecting things together (equations ( p8| ) and ( p0| ) ) we find that within distribu- 
tions 



1 



af'(x )(S + 5). 



(34) 



which gives us imediately the distributional limit of the regularized solution of the 
third equation of system (23), i.e., 



1 



N?(u) -» -af(x )(u + +6(u)). 



(35) 



Now we face the most complicated equation of system (|23|) which contains the 
distributionally ill-defined expressions, namely 

N:(u) = 2[N?(u)f'(x £ (u))p e (u)Y + [N?(u)f(x e (u)) Pe (u)}" 

- N?(u) f"{x t {u)) [9 * x t ] 2 {u) p e (u) - JVf («) f(x e (u)) p e (u) 
-N?(u)f'{x £ (u))x e (u)p £ (u), (36) 

where we have inserted all the dependences explicitly. Note that in addition to the 
obvious troublesome fifth term which involves the "square" of the ^-distribution 
(x e cx p e !), also the fourth term contains divergent and regularization-dependent 
contributions, since the derivatives of f{x) and N £ produce 5- and step functions 
respectively. We are going to calculate the distributional limits of the right hand 
side of equation ( [36] ) term by term, which can easily be done for the first three ones 
using yet well known techniques. 
First, second and third term: 

lim 2[N*{u) f'(x £ (u)) Pt (u)Y=laf(x ) 2 8 , 
lim [N?(u) f(x € (u)) p £ (u)}" = af(x )5 , 

c— >0 

lim N?(u)f"(x e (u))[6*x e } 2 (u)p e (u) =±af(xo) 2 f"(x )5 (37) 



We now come to the expressions involving the divergent and regularization- 
dependent factors. We outline the calculations in some detail. 
Fourth term: (iv) := —N®(u) f'(x e (u)) Pt {u) 
By integration by parts we obtain (tp again a test function) 

£ 1 

ip(u)N*(u)f'(x e (u))p e (u)du = [ y{eu)N?{eu)f{x e {eu))p{u)du (38) 



A 

1 1 



+ / tp(eu) Nf (eu) f" (x € (eu)) x € (eu) p(u) du + J tp(eu)N^(eu) f'(x e (eu))p(u)du. 

-l -l 



B C 



Inserting 9*9* Nf according to equations (|26j) and (|27]) for Nf we find for the 
expressions labelled A and B 
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A^-±af(x ) 2 S , (39) 
B^±af(x ) 2 f"(x )S. (40) 



Finally we have for the most complicated expression 

1 u 

C '= —a / ip(eu) f'(x t (eu)) p(u) / /' \x € (es)) p(s) ds du 



Ci 

l 

+ o° / <P( eu ) f'i x <i(eu)) plu) I f'(x e (es))sp(s)dsdu (41) 



2 

-l 



C2 



p(eu)/'(a: e (eu))p(u) / [a/'(a: e (e«)) + ijV?(ea)/"(i e (e«))]p(a)dadu 



C3 

The term labeled CI diverges due to the factor €~ x . However, as we shall see later, 
this expression combined with another one arising from the fifth term of equa- 
tion ( pq ) will provide us with a finite result. 

The term labelled C2 involves a regularization-dependent factor, i.e., 

1 u 

C2 ^af(x ) 2 S J p(u) J sp(s)dsdt. (42) 

-l -l 

As in the case of the previous term CI, this problem will be resolved later on, so 
that we are now left with the task of computing the limit of C3. Again using similar 
techniques we obtain 

Ci -> \ a/'(x ) 2 6 + ± af(xo) 2 f"(x ) S . (43) 
Collecting together the results of equations ([39]) , J40| ) and (^) , we finally have 

(iv) —C1 — C2 — > — i af'(x ) 2 5 + ~ a/'M 2 /"(io) <5 + 5 af'(x ) 2 S , (44) 

where now, by an abuse of notation CI and C2 denote the respective quantities in 
equation (|4l| ) without integration over the test function. 

Now we finally come to the last term in equation ( |36| ) which will resolve the problems 

with the expressions CI and C2. 

Fifth term: (v) := -N?(u) f'(x t (u)) x e (u) p t {u) 

Inserting solution (^5| ) we find 

<p(u)N?(u) f(x e (u))x e (u) p e (u) du 

1 1 

~Y a / f( eu ) f'i x e(eu)) 2 p(uf du p{eu) f'(x e (eu)) 2 up(u) 2 du . (45) 
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We see now that the expression (*) diverges, but that the sum CI + (*) converges. 
More precisely, by integration by parts we get 

1 us 

Cl + (*) = -\a J <p(eu)f'(x e (eu))p(u) J f"(x e (es))p(s) J f'(x £ (er))p(r)drdsdu 

-1 -1 -1 

1 U S 

^-i aV ?(0)/'(zo) 2 f"(xo) f P (u) f p(s) I p(r)drdsdu 



1/6 

= -^a<p(0)f(x ) 2 f"(xo). (46) 

Hence we have the distributional limit 

Cl + (*) - -±af'(xo) 2 f"(x )6. (47) 

On the other hand expression (**) contains a regularization-dependent factor which 
is proportional to the one contained in the term C2. More precisely, from equa- 
tion ( fl2| ) and the fact that 

l 

(**) -> -^af'(x ) 2 S J up{ufdu (48) 
-l 

we have 

1 u 

C2 + (**) 1 a f'(x ) 2 S [ f p{u) [ f sp{s) ds - up{u)] du] 



1/2 



= -\af(x ) 2 S. (49) 

We are now ready to collect together all five terms of equation ( p6|) , i.c^, of the 

second equation of the regularized system (p3|). From equations plj), (|44|), ( ^7| ) 
and ([l9]) we get the following result: 

N: -» a /(*„) 5 + j a /'(xq) 2 + (50) 

which, using initial condition (|24|), we readily integrate to 

JV e » -> 6 (! + «) + a [/(so) i/'(x o ) 2 (0H+ u+ )] . (51) 



Finally we have for the whole system (]23|) the following (distributional limits of the) 
solutions 

N?(u) = a(l + u) , 

N2(u) -^b(l + u) + a [f(x Q ) 5 + ^f'(xo) 2 (6(u) + u+) ] , 
N*{u)^) 1 af(x ){6(u)+u + ) , 

JV*(u) = 0. (52) 
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Hence, viewed distributionally, the Jacobi field suffers a kink, a jump and a (5-like 
pulse in the v-direction as well as a kink and jump in the cc-direction overlapping 
the linear flat space behavior. These effects can be understood heuristically from 
the corresponding behavior of the geodesies, given by equation (|l8|). The constant 
factor a, which gives the "scale" of all the nonlinear effects, arises from the "time 
advance" of the "nearby" geodesies, represented by the initial velocity of the Jacobi 
field in the it-direction (cf. initial conditions fl24|) ). 

Note, however, that this "time advance" is not the only effect generically generating 
deviations from the flat space behavior, but rather arises as an artifact of our initial 
conditions. One can easily show that different initial conditions on the deviation 
vector field, even without "time advance" , produce kinks and jumps as well. For 
example, let 

7V e a (-l) = (0,0,0,0) and N?(-l) = (0,0,6,0); (53) 

then a similar (but now even simpler) calculation leads to the following (distribu- 
tional limits of the) regularized Jacobi field 

iV £ " = , 

N v e - o(l + u) + \ b f'(x Q ) f"(x ) u+ + b f'(x ) 9{u) , 
Nf ^b(l + u) + hf"(x )u + , 

N* = . (54) 

The kink of the x-component of the deviation field now arises from the fact that 
a "nearby" geodesic passes the shock at an x-value of xq + b, hence according to 
equation ( |T8| ) suffers a kink of "strength" (1/2) f'(xo + 6). Taylor expansion yields 
f'(xo + b) « f'{xo) + f"{xo) b, such that the kink difference of "nearby" geodesies 
is given by (1/2) f"(xo) 6, which is exactly the factor given in the third equation 
of (|54]). The kink and jump in w-direction can be explained by similar heuristic 
arguments. 



IV. CONCLUSION 



We have shown that the geodesic and Jacobi equation of impulsive pp-waves can 
be treated consistently using the distributional form of the space time metric (||) . 
The main problem was to deal with products ill-defined within Schwartz's linear dis- 
tribution theory. The geodesic equation involves terms proportional to "65" while 
the Jacobi equation, although linear in the components of the deviation vector 
field, contains even more singular terms like "(5 2 " and "9 2 5". Our strategy con- 
sisted essentially in a careful and general regularization procedure, replacing the 
^-distribution by a mollifying sequence p e . While mathematically sound, this proce- 
dure corresponds to the physical idea of viewing the impulsive wave as the limiting 
case of a sandwich wave of ever decreasing support but constant strength. Our 
approach leads to smooth solutions of the regularized equations which possess a 
regularization-independent limit within the vector space of distributions. However, 
this is not obvious (cf. for exampleO) , but rather should strengthen our trust in 
impulsive waves as reasonable solutions of Einstein's equations. Furthermore our 
distributional "solutions" of the geodesic and Jacobi equation perfectly coincide 
with physical expectations, showing that the geometry of impulsive pp-waves can 
be described in a consistent way using the distributional form of the metric, as 
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long as one applies proper regularization methods as opposed to "ad-hoc strate- 
gies" involving certain "multiplication rules" . In fact when dealing with nonlinear 
operations (of a certain complexity) on singular, i.e., distributional data, reliable 
results can only be achieved by careful regularization procedures (cf. thexomments 
on "common errors" even in the mathematical literature given by HajeklB). 
However from an even more formal point of view, the approach taken in this pa- 
per suffers from the fact that we do not really have an elaborated solution concept 
for the regularized equations ll ^] ir '^ a notion is provided by the theory of Nonlin- 
ear Generalized FunctionsEj^Lj due to J. F. Colombeau and others, where - loosely 
speaking - the regularized sequence is viewed as one mathematical object, i.e., a gen- 
eralized function. Recently, singular, nonlinear ODEs (exactly the type of equations 
we were deaJiig with) have been studied in the context of this new mathematical 
fraineworkt3'E3, providing a powerful tool to handle both, singular data and espe- 
cially this type of nonlinear differential equations. Future work will be concerned 
with the analysis of the geodesic and Jacobi equation for pp-waves from the point 
of view of this new solution concept. 
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APPENDIX: 



In this appendix we prove (for small e) the global existence and uniform boundedness 
(in e) of the solutions to the last equation of system ([?]). For simplicity we only treat 
the one-dimensional case and write g instead of Df. The main idea is to use || 
bounds on the regularizing sequence p t . We start with the following statement. 

Proposition: Let g : K — > R smooth and p € T>([— 1, 1]), J p = 1, p e (x) :— 
(l/e)p(x/e). For fixed e consider the differential equation 

x e (t) = g{x e (t)) Pe {t) , (Al) 

with initial conditions x e (—l) — xo and x £ (— 1) = Xq > 0. Then there exists a 
unique solution x e (t) on the interval J e := [— 1, — e + a], where 

ttt 1 (b > 0) , (A2) 



Mloo,/ \\p\\i + x ' 2L 
I := {x £ R : | x — xq \< b + xq } and L is a Lipschitz constant for g on I . 

The Proof is just a modification of the proof of the classical (first order) ODE 
existence and uniqueness theorem. We work on the (nonempty, closed) subset X e := 
{x e S C(J C ) : | x e (t) — x \< b + x } of the Banach space of continuous functions 
on the int erva l J £ . Since no first order derivative of x e enters the right hand side of 
equation (|Al|) we may define the integration operator A by 

t s 

[Ax e ](t) := x + x (t + l)+ / / gix^r)) p e (r) dr ds . (A3) 
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By construction Ax t G C(J e ), but we even have Ax e G X e since 

t s 

I [Ax e ]{t)-x | <x (t + l)+ [ [ \g(x t (r))\\ Pt (r)\ drds 



<a(x + \\g\\oo,i\\p\\i)+x <b + x (A4) 

By a similar estimate one shows that A, in fact, is a contraction on X e . Hence by 
the fixed point theorem we have a unique solution x e G X e . □ 

Since a is independent of e we have f or s mall e: — e + a > e. But at t = e the right 



hand side of the differential equation (Al) has already been "turned off" . Hence for 



small e the solutions are defined on the whole real line. 

Furthermore, since for all e the solution lies in X e , x e is bounded uniformly in e 
for say t < a/2 (for e small enough). For larger values of t x t grows only linearly 
with x e which is bounded by io + \\g\\oc.i\\p\\i- Hence (for small e) x e (t) is bounded 
uniformly in e on every compact set. 
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